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ABSTRACT _ 

We develop and apply the position angle-only shear estimator of |Whittaker et al.| ( |2014| 
to realistic galaxy images. This is done by demonstrating the method on the simulations 
of the GREAT3 challenge (Mandelbaum et al. 2014b| ), which include contributions from 
anisotropic PSFs. We measure the position angles of the galaxies using three distinct methods 
- the integrated light method, quadrupole moments of surface brightness, and using model- 
based ellipticity measurements provided by IM3SHAPE. A weighting scheme is adopted 
to address biases in the position angle measurements which arise in the presence of an 
anisotropic PSF. Biases on the shear estimates, due to measurement errors on the posi¬ 
tion angles and correlations between the measurement errors and the true position angles, 
are corrected for using simulated galaxy images and an iterative procedure. The proper¬ 
ties of the simulations are estimated using the deep field images provided as part of the 
challenge. A method is developed to match the distributions of galaxy fluxes and half- 
light radii from the deep fields to the corresponding distributions in the field of interest. 
We recover angle-only shear estimates with a performance close to current well-established 
model and moments-based methods for all three angle measurement techniques. The Q- 
values for all three methods are found to be Q ~ 400. The code is freely available online 
at http://www.jb.man.ac.uk/-mbrown/angle_only_shear/ 

Key words: gravitational lensing: weak - methods: analytical - methods: statistical - cosmol¬ 
ogy: theory 


1 INTRODUCTION 

Gravitational lensing provides a way of probing the dark Universe. 
As light travels to us from distant objects, it is deflected by fore¬ 
ground mass inhomogeneities. This effect can be exploited to pro¬ 
vide information on the dishibution of dark matter and the nature 
of dark energy, enhancing our understanding of the late-time evo¬ 
lution of the Universe. In the weak regime, the lensing of light by 
the large-scale structure of the Universe induces small but coherent 
distortions of observed background galaxy shapes. Since one does 
not know the intrinsic galaxy shapes, statistical methods must be 
employed to extract information about the underlying mass distri¬ 
bution. 

If we assume that the size of the source object is much smaller 
than the angular scale over which the gravitational potential of the 
foreground object changes, the distortion of the image can be de¬ 
scribed by a linearized mapping of the object from the source plane 
to the image plane (see [Bartelmann & Schneider] l |200l| > for an in- 
depth discussion). This linear transformation can be separated into 
two components; the convergence, which magnifies the image, and 
the shear, which both magnifies and changes the shape of the im¬ 
age. In this paper, we are concerned with estimating the shear from 


a field of simulated galaxy images. Assuming that the shape of a 
source galaxy can be described by an ellipse, the shear changes the 
ellipticity of the galaxy in the image plane. If one assumes that the 
distribution of the intrinsic ellipticities (i.e. the ellipticities in the 
source plane) are randomly orientated, then an accurate measure¬ 
ment of an individual galaxy’s ellipticity provides a noisy, unbi¬ 
ased estimate of the shear at the position of the galaxy. This is the 
standard method for performing a weak lensing analysis. Such an 
approach has already been used to constrain various cosmological 
parameters, such as the amplitude of the matter power spectrum 


(e.g. |Brown et al.|2003j |Hoekstra et al.|2006[ Fu et al. 2008 


the dark energy equation of state (e.g.fSchrabback et al. 2010 |Kil-| 

|binger et al. |20 1 3) . 


and 


In order to perform weak lensing using the standard method, 
one requires accurate measurements of the galaxy shapes. There are 
generally two classes of shape measurement techniques i j Viola et ah| 
|2014[ >; these are the moments-based and the model-based methods. 
Moments-based methods measure the ellipticity of a galaxy by cal¬ 
culating the quadrupole moments of the light distribution of the 
pixelized galaxy image. Examples of such algorithms include the 
Kaiser-Squires-Broadhurst (KSB) metho d j Kaiser et al.|1995| l and 
the DEIMOS method ^Melchior et al.pOll) . Model-based meth- 
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ods fit a parameterized galaxy model to a pixelized galaxy image 
by finding the extremum of a loss function. Examples of this ap¬ 
proach include lensfit l Miller et al.|2007[|Kitching et al.|2008) 
and IM3SHAPE l |Zuntz et al.|2013| >~ 


For any realistic weak lensing survey, measurements of the 
galaxy shapes need to be recovered from noisy, pixelized images. 
These images are convolved with an instrumental and/or atmo¬ 
spheric point spread function (PSF), and both moments and model- 
based methods must correct for the PSF contribution. Moments 
based methods generally also implement a weighting function to 
reduce the effects of noise at large scales, and this function must 
subsequently be corrected for. For the case of an isotropic PSF, an 
incorrect calibration of the required correction will lead to multi¬ 
plicative biases in the shear estimates. However, the orientation of 
a galaxy will remain unaltered. The same is true if one multiplies 
the galaxy image by a circular weighting function centered on the 
centroid of the galaxy. In |Whittaker et al.| ( (2014[ > (hereafter W14), 
we proposed a method of performing weak lensing using only mea¬ 
surements of the galaxy orientations with the aim of exploiting this 
property. We demonstrated the method using simple simulations 
where we ignored the effects of a PSF and considered a simple 
Sersic galaxy model for the intrinsic galaxy shapes. We found that 
the position angle-only method has the potential to yield shear es¬ 
timates with a performance comparable with the KSB method. 

An anisotropic PSF will bias position angle measurements 
and, therefore, estimates of the shear if not corrected for. This bias 
is addressed in this paper where we build upon the ideas introduced 
in W14 and apply the angle-only method to the simulations of the 
control-ground-constant (cgc) branch of the GREAT3 simulations 
l |Mandelbaum et ah||2014b| . These simulations were designed to 
test the performance of a shear estimator on realistic galaxy images 
and include the effects of noise, pixelization, an anisotropic PSF, 
and realistic distributions of galaxy flux, size and shape. 

In Section [2] we summarize the angle-only method proposed 
in W14 and discuss how noise on the position angle measurements 
impacts the analysis. We describe the three methods used to mea¬ 
sure the position angles of the galaxies in Section[3] The application 
of the angle-only method to the GREAT3 simulations is presented 
in Section]?] where we describe the weighting scheme used to mit¬ 
igate the effects of an anisotropic PSF and the calibration simula¬ 
tions implemented to remove the effects of the biasing introduced 
by measurement errors on the position angles. The results are pre¬ 
sented in Section]?] where we compare the angle-only method, us¬ 
ing the three angle measurement techniques, with the results of a 
naive application of IM3SHAPE and also with the highest entries 
to the GREAT3 challenge from IM3SHAPE and the KSB method. 
We conclude with a discussion in Section[6] 


2 THE ANGLE-ONLY METHOD 


In this section, we briefly describe the angle-only shear estimator 
introduced in W14. 

Working within the regime of weak lensing, we can express 
the observed complex ellipticity of a galaxy, e obs , in terms of the 
intrinsic ellipticity, e lnt , and the reduced shear, g: 


^int 1 „ 

^obs _ e + 9 
1 + g*e lnt ’ 


(1) 


where the asterisk denotes complex conjugation. We can also ex¬ 


press e obs in polar form, such that 

e° bs = | e obs | cos (2a), 

e° bs = | e obs | sin (2a), (2) 

where a is the observed position angle of the galaxy. If we as¬ 
sume a constant shear, and that the intrinsic ellipticities are drawn 
randomly from an azimuthally symmetric probability distribution, 
/ (e mt ). the means of the cosines and sines of the position angles 
can be written as 


(cos (2a)) =Fi (|g|) cos (2a 0 ), 

(sin (2a)) =Ei (|g|) sin (2a 0 ), (3) 


where ao is the shear position angle. The function Ei (|g|) is found 
to be 


(|fl|) = l ^ K ” X| f ] da int d |e int | / (|e int |) 

x hi (|g|, |e lnt | ,a lnt ) , (4) 

where a lnt denotes the intrinsic position angle. The function 
hi (|g| , |e lnt | , a mt ) is given as 


with 

Ei 

E 2 


1 (I I I int I int \ Eh 

Isl (l + |e in f) + (1 + Is! 2 ) |e int | cos (2a int ) , 
(l — ls| 2 ) | e ' nt | sin (2a int ) . 


(5) 


( 6 ) 


From here, W14 showed that, for a finite number of source galaxies, 
an estimate of the orientation of the shear is given by 



/ Ef =1 sin(2qW) \ 

VEf = iCOs(2aW)J ’ 


(7) 


where N is the number of galaxies in the sample. A preferred value 
of the Fi (|g|) function can also be recovered for the sample using 



r n 1 

2 

r 1 n 1 

Ei (Isl) = ^ 

— ^ c °s(2a (i) ) 
i=l 

+ 

— ^sin( 2 a (i) ) 
i= 1 


(8) 

This relation can then be inverted to provide an estimate of |g|. 
We can combine this estimate with an estimate of do provided by 
equation 0 to recover a complete estimate of the shear. 

Measurement errors on the position angles bias estimates of 
the mean trigonometric functions in equations 0 and 0- This bias 
propagates into estimates of the shear. As discussed in W14, for a 
general error distribution we can write the measured position angle 
as 


d = a + 5a, (9) 

where 5a is the error on the measurement. The means of the 
trigonometric functions are then 

(cos (2d)) = (cos (2a + 25a)), 

(sin (2d)) = (sin (2a + 25a)) . (10) 
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We can write equation ( | 1 0| ) in terms of the covariance between the 
true position angles and the measurement errors, such that 

(cos (2d)) =C' + cov (cos (2a), cos (25a)) 

— cov (sin (2a), sin (25a)), 

(sin (2d)) —S' + cov (sin (2a) , cos (25a)) 

+ cov (cos (2a), sin (25a)), (11) 

where we define 

C' = (cos (2a)) /3 C — (sin (2a)) /3 S , 

S' = (sin (2a)) /3 C + (cos (2a)) /3 S , 

/ 3 C = (cos (25a)), 

[} s = (sin (25a)). (12) 

Hence, we see that errors on the position angles of the galaxies and 
correlations between the true lensed position angles and the mea¬ 
surement errors bias our estimates of the mean trigonometric func¬ 
tions. This in turn biases the shear estimates. W14 addressed this 
issue, in the absence of a PSF, by estimating the f3 and covariance 
terms using simulations and applying an iterative procedure. How¬ 
ever, we find that this approach is inadequate for fields which have 
a large anisotropic PSF. In Section[4] we introduce an alternative it¬ 
erative procedure. We use an initial estimate of the shear recovered 
using weighted averages of the trigonometric functions to simulate 
the expected bias in the shear estimates and subsequently correct 
our estimates. 


sky. This step is iterated until the difference between subsequent 
estimates of the components of the centroid are less than 1CP 3 of a 
pixel. 

We then estimate the ID integrated light distribution, /' (9), of 
the galaxy as a function of assumed position angle, 9, by integrating 
I' (0) over the radial component centered on the estimated centroid. 


I'(9) = J drI w (r, 9). (14) 

The details of how this integral is performed are discussed in W14. 
Finally, the estimated position angle of the galaxy is given by 


„ _ 1 , _i ffdOr (0) sin (20) ^ 
a ~ 2 tan V /dW' (9) cos (29)) ' 


(15) 


3.2 A moments based method 


The moments based method uses second order moments of the con¬ 
volved weighted brightness distribution, I w (0), to estimate the po¬ 
sition angle ( jKaiser et al.|1995| . We follow the approach outlined in 
the previous subsection to estimate the width of the weighting func¬ 
tion used for each field and the centroid of the individual galaxy im¬ 
ages. The position angles of the galaxies are then calculated using 


the second-order moments of (0), as defined in equation i 13 i 
such that 


.. 1 4 -if 2Qu 

a = - tan —-—— 

2 \ Q 20 — Q 02 


(16) 


3 MEASURING THE POSITION ANGLES 

When analysing the cgc branch of the GREAT3 simulations, we 
use three distinct methods to measure the position angles from the 
simulated galaxy images. The first is the integrated light method 
discussed in W14. The second uses the second order moments of 
brightness of the galaxy image. The third uses ellipticities mea¬ 
sured by IM3SHAPE to determine the position angles. Here we 
briefly discuss the three approaches. 


3.1 The integrated light method 


The integrated light method for measuring galaxy position angles 
is discussed in detail in W14. Here we discuss the application of 
the method to the images of the cgc simulations. 

We begin by estimating the mean half-light radius of the 
galaxies in the field of interest using the relationship between the 
mean half-light radius and the mean flux of the galaxies discussed 
in Section[4]and given in equation ( |2 1 1 ). For each galaxy in the cgc 
branch, we initially assume the centroid of the galaxy to be the 
centre of the image. We apply a circular Gaussian weighting func¬ 
tion to the image centered on this initial estimate of the centroid. 
The half-light radius of the weighting function is equal to twice 
the mean half-light radius of the galaxies in the field. The image is 
then convolved with a circular Gaussian kernel with a width of two 
pixels to reduce the effects of pixelization. The centroid is then re¬ 
calculated using the first-order moments of the convolved weighted 
surface brightness distribution, I w (0), with the moments defined as 


Qij. 


..k 


f d 2 9I w (0)9i9j...9k 
Jd 2 9I w (0) 


(13) 


and where 0 is the angular position of the galaxy image on the 
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3.3 Using I M3 SHAPE 


Using IM3SHAPE, we measure the observed ellipticities of the 
galaxies in each held and use the ellipticity measurements to de¬ 
termine the galaxy position angles. When fitting a model of the 
galaxy shape to the galaxy image, IM3SHAPE takes into account 
the contribution of the PSF by convolving the model galaxy with 
a star held image. The GREAT3 challenge provides star held im¬ 
ages at the same resolution as the galaxy images for each held of 
the cgc branch. When using IM3SHAPE to measure the position 
angles of the galaxies, we use these star held images directly to 
calibrate for the PSF. However, when we require full ellipticity in¬ 
formation (i.e. when estimating / (|e lnt |) in Section [ 4 ] and when 
using IM3SHAPE to estimate the shear, as a comparison with the 
angle-only method, in Section [5}, we use star held images which 
are upsampled by a factor of seven using Gal Sim ( |Rowe et al.| 
|2014[ >. We use original star held images for the position angle mea¬ 
surements to reduce computation time when calibrating the angle- 
only shear estimates using simulations. Since IM3SHAPE corrects 
for the PSF when measuring the ellipticities, the angle-only method 
using IM3SHAPE differs from the integrated light and moments 
based methods which do not include any PSF correction at this 
stage. Upon recovering measurements of e obs , we estimate the po¬ 
sition angles of the galaxies as 




(17) 


4 APPLICATION TO THE GREAT3 SIMULATIONS 

The cgc branch of the GREAT3 challenge simulates ground-based 
observations of 200 10x10 deg 2 helds. Each simulated observation 
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contains 10 4 resolved galaxy images. A constant shear is applied to 
all of the galaxies within a particular field. The images of the galax¬ 
ies within each field are convolved with a constant anisotropic PSF, 
and a constant level of background Gaussian noise is assumed. The 
applied PSF and background noise levels are varied between the 
different observations, and the applied shear is different for each 
field observed. Noiseless star field images are included for each 
field providing an image of the PSF. The cgc branch also includes 
five deep field observations for use as a training dataset. These ob¬ 
servations consist of galaxy images that are one magnitude deeper 
than the challenge observations, but the dataset retains only images 
of the galaxies which would be present in the rest of the challenge. 
The properties of the simulations are discussed in detail in |MandeT| 
|baum et aL] < |2014b) . 

In this section, we describe the application of the angle-only 
method to the cgc branch of the GREAT3 simulations. We mea¬ 
sure the position angles of the galaxies from the simulated images 
using the three methods discussed in the previous section. We begin 
by outlining the steps of the procedure followed in the angle-only 
analysis. 

(i) Calculate the F\ (|g|) function using equation 0. This re¬ 
quires an estimate of / (| e lnt |) which we obtain from the GREAT3 
deep field calibration sets using IM3SHAPE. 

(ii) Measure the fluxes of the galaxies in the deep field images 
by summing over the pixel values and recover measurements of the 
half-light radii and the bulge to total flux ratio of the galaxies from 
IM3SHAPE, to be used in calibration simulations. 

(iii) Determine the relationship between the mean fluxes and the 
mean half-light radii of the five deep fields. 

(iv) Measure the mean flux of the galaxies in the field of interest 
and use this measurement to modify the fluxes and half-light radii 
of the deep field measurements output from step (ii) above. 

(v) Use the modified fluxes and half-light radii with the intrin¬ 
sic ellipticity estimates, bulge to total flux ratio measurements, and 
a suite of random uniform distributed intrinsic position angles to 
simulate zero shear galaxy images using GalSim. The noise in 
the simulated images is estimated from the GREAT3 image being 
analysed. 

(vi) Measure the position angles of the galaxies in the zero shear 
simulations using one of the three methods discussed in the previ¬ 
ous section. These measurements are used to construct a weighting 
function the purpose of which is to correct for PSF anisotropy and 
pixelization effects. 

(vii) Measure the position angles of the galaxies in the field 
of interest and estimate the shear using weighted averages of the 
trigonometric functions. 

(viii) Use the shear estimates to produce an updated set of sim¬ 
ulations which include information about the shear. 

(ix) Estimate the shear in the updated simulations using the 
same weighting function as in step (vi). Use these estimates to de¬ 
termine the bias in the estimates recovered in step (vii) and correct 
for the bias. 

(x) Repeat steps (vii)-(ix) using the corrected shear estimates as 
the input shear for the simulations until the estimated bias is below 
a desired threshold value. 

We now discuss our application of this procedure to the cgc 
branch in detail. 

We began by estimating the form of the F\ (|gj) function (step 
(i)) using IM3SHAPE to measure the observed ellipticities of the 



Figure 2. The distribution of position angles measured from field 2 using 
the integrated light method (black curve). The vertical black line shows the 
position angle of the PSF. Here we see that the measured position angles 
are biased in the direction of the PSF. The red curve shows the distribution 
of measured position angles recovered from the zero shear simulations. We 
see that the distribution from the simulations provides a good description of 
the distribution from the image data, with the difference between the two 
being attributed to the shear signal for the zeroth-order shear estimate. 


galaxies for each of the deep fields of the cgc branch. When mea¬ 
suring a galaxy’s ellipticity using IM3SHAPE, we used star im¬ 
ages which were upsampled by a factor of seven to correct for the 
PSF contribution (as opposed to the original star images used when 
we were concerned only with a galaxy’s position angle). For each 
of the deep fields, we averaged over the observed ellipticity mea¬ 
surements to estimate the shear for that field. Using the estimated 
shears and the measured ellipticities, we recovered estimates of the 
intrinsic ellipticities of the galaxies by inverting equation 0 and 
reconstructed / (|e lnt |), shown in Figure 1 From this estimate of 
the distribution, we calculated the F\ (|g|) function numerically 
using equation 0; *is is shown in the right panel of Figure |T| 
As discussed in W14, the integral carried out when calculating the 
F\ (|g|) function smooths the / (|e lnt |) distribution. However, if 
we use a large binsize for the distribution we lose information. We 
therefore chose a binsize of 5 x 10~ 3 in accordance with the binsize 
used in W14. 

The effects of the PSF on a galaxy’s observed ellipticity and 
orientation depend on the underlying ellipticity of the galaxy. When 
using the angle-only estimator, we do not recover information about 
the ellipticities, and therefore a complete understanding of the ef¬ 
fects of the PSF on an individual galaxy’s orientation cannot be 
realised. Instead, we chose to examine the impact of the PSF on 
the ensemble of galaxies. In the absence of a PSF (and pixeliza¬ 
tion), we expect the distribution of position angle measurements to 
be approximately uniform, with any deviation from uniformity be¬ 
ing attributed to the underlying shear signal. In Figure [2] we show 
the distribution of measured position angles from field 2 of the 
GREAT3 cgc simulation set measured using the integrated light 
method. We see the unsurprising result that the position angles are 
biased in the direction of the PSF. This results in non-zero corre- 
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Figure 1. Left-hand panel: The reconstructed f (|e int |) using IM3SHAPE to recover the estimates of the intrinsic ellipticities of the 
galaxies in the deep field images of the cgc branch and with a bin width of A |e int | = 5 X 10 —3 . Right-hand panel: The F\ (|gr|) 
function calculated using equation jJj and the reconstructed / (|e int |) shown in the left-hand panel. 


lations between the galaxy position angles and the errors on the 
measurements which biases the shear estimates via equation O' 
To correct for this effect, we adopt a weighting scheme in 
which we downweight the contribution to the mean trigonomet¬ 
ric functions (see equations 0 and {8}) from galaxies which align 
with the PSF. To understand how the galaxies should be weighted, 
let us assume a large sample of galaxies with a zero shear signal. In 
the absence of a PSF, we expect a uniform distribution of measured 
position angles and therefore the mean trigonometric functions to 
be zero. The observed unit vector of the galaxy is 

( cos (2a) \ 

n ^=Un(2d) )■ (18) 

In the presence of a PSF, the angle distribution becomes non- 
uniform, and the mean unit vector, (n), will be 

(n) = f dd n (a) /psf (d), (19) 

“'“f 

where /psf (d) is the distribution of measured position angles 
given a non-zero PSF. If we introduce a weighting function, w (d), 
such that 

( n } = / dd w (d) n (d) /psf (d) , (20) 

it is clear that one can correct for the effects of the PSF, such that 
we recover mean trigonometric functions equal to zero, if we set 
w (d) = 1 //psf (d). To proceed, we therefore require an estimate 
of the distribution of measured position angles when a PSF is in¬ 
cluded. 

To achieve this, we used GalSim to simulate 10 5 galaxy im¬ 
ages with specific properties provided by IM3SHAPE and assum¬ 
ing a zero input shear signal. For each galaxy in the five deep fields, 
we measured e° , the bulge to total flux ratio ( B/S ), and the half- 
light radius (R e ). For each field, we also estimated the shear by 
averaging over e obs allowing us to estimate the e lnt of each galaxy 



<R e > (arcsec) 


Figure 4. The mean flux of the galaxies plotted as a function of the mean 
half-light radius for each of the five deep fields. Here we see a linear rela¬ 
tionship between the two quantities. The line is the best fit to the data. 

by inverting equation 0 - For each of these galaxies, we simulate a 
corresponding galaxy using the measured properties as a template. 
The intensity profile of each simulated galaxy was assumed to be a 
Sersic profile consisting of a bulge with Sersic index n s = 4 and 
a disc with n s = 1, and with the half-light radii of each of these 
components being identical. The simulated galaxy is created with 
an intrinsic ellipticity of | e mt | and with a position angle drawn ran¬ 
domly from a uniform distribution with the range [—90°, 90°). 

For each field in the cgc branch, the galaxy images are con¬ 
volved with a different PSF, and only resolved galaxy images are 
included in the sample. Including only resolved galaxy images im- 
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Figure 3. The top left-hand panel shows the distributions of galaxy fluxes for the five deep fields of the cgc branch. The top right-hand 
panel shows the distributions of half-light radii for the galaxies in these fields, with the half-light radii measured using IM3SHAPE. In 
the bottom two panels we have modified the distributions of fields 1-4 using the relationships given by equations <22|25t so that they 
match the distributions of field 0. These plots indicate that the flux and half-light radii distributions in a particular field can be accurately 
reproduced using these relationships. 


plies that the distributions of the half-light radii and, therefore, 
the fluxes of the galaxies are different for each field. In order to 
construct accurate calibration simulations, we require the distri¬ 
butions of these properties in the simulations to match the distri¬ 
butions in the field being analyzed. We used the properties esti¬ 
mated by I M3 SHAPE when analysing the deep field images to in¬ 
vestigate how this can be achieved. First, we measured the flux of 
each galaxy by summing over the pixel values in each galaxy im¬ 
age, and we recovered the half-light radii estimates provided by 
IM3SHAPE (step (ii)). Figure [3] shows the distributions of galaxy 
fluxes and half-light radii for each of the five deep fields. We see 
the difference between the distributions for each field. In Figure 
[4] we plot the mean galaxy flux for each field as a function of 
the mean half-light radius. This indicates a linear relationship be¬ 


tween the two quantities. We fitted a linear function of the form 
(S) = m (R e ) + c to this data (step (iii)) with the fitted parame¬ 
ters given by m = 162.4 ADU arcsec -1 and c = —14.19 ADU. 
This relationship allows for a direct estimate of (R e ) simply by 
summing over the pixel values in a particular field, such that 

(R e ) = (21) 

m 

One should not assign significance to the relationship shown in Fig¬ 
ure [ 4 ] The method for recovering estimates of ( R e ) works as long 
as one can characterize a relation between (R e ) and (S). 

If Pg ( S ) is the distribution of the fluxes in field n and (S) 
is the mean flux of a galaxy for that field, we find that we can write 
an approximate relationship for the flux distributions of two fields 
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a and b as 


that 


Pg(S) = P b s (KS), (22) 


where 


K = 


(SL 

(s ) b ■ 


(23) 


We find we can also write a similar relationship for the distributions 
of half-light radii of the two fields, ( R e ) and ( R e ), 


Pg 0 (P o ) = Pl(Pe + A W ), (24) 


where 


A" = (R e ) a - (R e ) b . (25) 

If one measures the means of the fluxes of the galaxies in fields 
a and b, one can recover estimates of the means of the half-light 
radii in these fields using the linear relationship shown in Figure [4] 
and equation ( |21fr . It is then possible to modify the distributions of 
the fluxes and half-light radii in field b so that they approximately 
match those in field a. An example of this procedure is shown in 
the bottom two panels of Figure [3] where we have modified the 
distributions of the deep fields 1-4 to match deep field 0. We see 
that the distributions of the fluxes and half-light radii of a particular 
field can be accurately reproduced if the mean of the galaxy fluxes 
can be accurately measured and, as this requires simply summing 
over the image pixels within a particular field, we expect this to be 
achievable. 

When applying this method to the challenge data, we mea¬ 
sured the mean flux of the galaxies in the field of interest. We then 
modified the fluxes and half-light radii of all the galaxies in the five 
deep fields. We combined the modified fluxes and half-light radii 
with the estimated values of |e mt | and B/S, and a suite of uni¬ 
form random intrinsic position angles to provide us with the inputs 
required for 5 x 10 4 simulated galaxies (step (iv)). The GREAT3 
simulations use galaxy pairs whereby for each galaxy with an in¬ 
trinsic orientation of a mt there is an identical galaxy with an in¬ 
trinsic orientation of a lnt + 90°. This is done to reduce intrinsic 
shape noise in the shear estimates and hence the number of galaxies 
needed to average over. For each of our 5 x 10 4 simulated galax¬ 
ies, we therefore created an identical galaxy image but rotated by 
90° with respect to the first, mimicking the procedure implemented 
by GREAT3. This provided us with 10° simulated galaxy images. 
These images were convolved with the upsampled star field images 
to simulate the effects of the PSF. The pixel scale of the simulations 
was chosen to match the pixel scale used by GREAT3. The noise in 
the images was assumed to be Gaussian with the variance estimated 
from the values of the outermost pixels of each 48 x 48 pixel 4 stamp 
in the field and with all 10 4 stamps used in the estimate (step(v)). 

Once the simulated galaxy images were created, we measured 
the position angles of the galaxies using each of the three meth¬ 
ods discussed above (step (vi)). The results of this for field 2 when 
using the integrated light method are presented as the red curve in 
Figure[2] This distribution was binned with a bin size of 1° to give 
the distribution of position angles, P (at). As described above, the 
required weighting function is the reciprocal of this distribution 

w{a ' ) = ph~y (26) 

We used this weighting function to correct the averages of the ob¬ 
served trigonometric functions for the PSF and pixelization, such 


(n) 


ZiLivfr) 


(27) 


Ideally, the weighting applied to each galaxy in the image would 
depend on the position angle of the galaxy without the effect of 
lensing (but with the effects of PSF, pixelization and noise still 
present). However, this is obviously not possible for real data. As¬ 
suming that the rotation induced by lensing is small, the weighting 
of the trigonometric functions will not be significantly affected by 
lensing and, provided the simulations are accurate, any resulting 
bias to the shear estimates should be corrected for by the iterative 
procedure discussed below. 

Using the weighted averages given by equation \21\ , we re¬ 
cover a zeroth-order estimate of the input shear signal. g <0 \ in 
each field via equations (J 7 J and |8| (step (vii)). Let m' and c[ re¬ 
spectively denote the multiplicative and additive biases on the com¬ 
ponents of the zeroth-order shear estimates (Heymans et ah 2006[ 
|Huterer et al.|2006||Massey et al.|2007) . The estimates can then be 
written as 


g? ] = Si + m'i9i +c'i + 8gf ] , (28) 

where 8g f 1 ' is an error on the estimate with zero mean. If we as¬ 
sume that the weighting function successfully mitigates the effects 
of an additive bias, the zeroth-order estimate is expected to underes¬ 
timate the modulus of the true shear signal as there is no correction 
for the measurement error bias arising from /3 C in equation ( |12| >. 
It is possible to correct for this bias using the zero shear simula¬ 
tions to model the measurement error for each field as discussed in 
W14. However, for cases where there are large contributions from 
anisotropic PSFs, we find that the errors on the position angle mea¬ 
surements can be large leading to small values of /3 C . Hence, at¬ 
tempting to correct for these biases can lead to substantial outliers 
in the shear estimates. Also, errors on the estimates of /3 C propa¬ 
gate nonlinearly into estimates of the shear. We, therefore, choose 
to simulate a further suite of 10 5 galaxy images using the procedure 
described above, and we shear the galaxies using the estimates g (0) 
(step (viii)) for each field. We measure the shear from these updated 
simulations and use these estimates to determine, and correct for, 
the bias in the initial estimates recovered for each field (step (ix)), 
such that 

9 W = 9 W - (siim - 9 (0) ) > (29) 

where the shear estimate recovered from the first-order simulations 
is 


s (i) _ a(0 ) 


9i 0) + m i9i°^ + G + 8g^ ] . 


(30) 


This method, therefore, corrects for the bias introduced to the esti¬ 
mates of the shear and does not correct for the bias on the individ¬ 
ual galaxy position angle estimates. We can iterate this step until 
estimates of the shear between subsequent iterations are consistent 
(step (x)). 

For a large number of simulated galaxy images, noise in the 
estimated bias will be sub-dominant to noise in the zeroth-order 
shear estimate. However, using this iterative method, noise in the 
simulated shear estimates, <^7 , due to a finite number of simu¬ 
lated galaxy images, propagates linearly into the final shear esti¬ 
mates. Therefore, if we assume that the details of the simulations 
are accurate, there is no additional noise bias expected from this 
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procedure. For the n th iteration, we can write 

9 (n) =9 (0) -(g^-9 (n - 1) )- (3D 

Assuming that the simulations provide an accurate description of 
the true field, we show in Appendix A that the residual bias for the 
n th iteration is 



which converges to zero for |m;| < 1. As explained above, the 
zeroth-order shear estimates are expected to underestimate the 
modulus of the true shear. Hence, we expect m' to be confined 
to the range — 1 < m'i < 0, and therefore the shear estimates to 
converge for all fields. 


5 RESULTS 


Here we compare the results of the angle-only approach using each 
of the three methods to measure the position angles of the galaxies 
discussed in Section[3] We include the results obtained from a naive 
application of IM3 SHAPE where full ellipticity information is used 
(measured using the upsampled star images to correct for the PSF) 
with no additional calibration scheme. 

The Q-value is the metric used to quantify the performance 
of the estimators in the GREAT3 challenge ( |Mandelbaum et ak| 
|2014b| l. It is defined as a function of the multiplicative and additive 
biases (m and c respectively ) where the biases satisfy the approx¬ 
imation (Heymans et al. 2006; Huterer et al. 2006; Massey et al. 
|2007l > 

gi = (1 + rrii) gi + a + Sgt. (33) 


The subscript i denotes the components of the shear in a reference 
frame aligned with the PSF in the field being analysed. The biases 
are estimated using a linear regression over the 200 fields given the 
known input shears. The Q value is then calculated as 


Q 



2000 ?) 


( mi ) +V._ ( Ci ^ 

y ^target J —+,X Retarget J 


(34) 


where the subscript + corresponds to the components of the biases 
in the direction aligned with the PSF for each field and the subscript 
x denotes the direction perpendicular to the PSF. The target values 
of the biases, m ta r g et and Ctarget, are based on the requirements of 
the ESA Euclid space mission {Massey et al.|2013) and are given 
as m t arget = 2 x 10 _,i and c ta rget = 2 x 10~ 4 . The constant r] 
normalizes the metric such that a value of Q ss 1000 is expected 
for estimates of the shear which achieve the target values of m and 
c. The term <7 2 in corresponds to the typical dispersion of the biases 
due to pixel noise and is determined using trial submissions to the 
GREAT3 challenge. For the challenge, the values adopted for the 
ege branch were r] = 1.232 and af nul = 4. We estimated the Q- 
values for our methods using the publicly available GREAT3 metric 
evaluation scripQ This script also provides the estimates of m, and 
d used to calculate the Q-value and their corresponding error bars. 


In the discussion that follows, we include the highest Q- 
value entry to the ege branch of the GREAT3 challeng^] from 
IM3SHAPE - for a comparison of the performance of the angle- 
only method with the challenge submission of a model based 
method. This entry implements a multiplicative calibration factor to 
correct for noise biases expected to arise in Maximum-Likelihood 
shape estimation. Additive biases are expected in the presence of 
anisotropic PSFs. However, no calibration for this bias is included. 
We also include the highest entry submitted to the challenge us¬ 
ing the KSB method - for a comparison with a moments based 
method. The details of the submissions using IM3SHAPE and the 
KSB method are presented in ( Mandelbaum et al.|2014a) . For the 
angle-only analyses, the results presented made use of a single it¬ 
eration of the procedure discussed in Section]?] 

Figure[5]shows the recovered angle-only shear estimates when 
using the integrated light and moments based methods to measure 
the position angles. These plots display the estimates for the com¬ 
ponents of the shear in a reference frame aligned with the simulated 
images. It should be emphasized that this is not the reference frame 
used to measure the Q-value, which employs a coordinate system 
that is aligned with the PSF for that field as described above. The 
error bars have been estimated using a linear approximation of the 
Fi (|g|) function as discussed in Appendix B. Figure [6] shows the 
results of I M3 SHAPE using both the angle-only method and full 
ellipticity information - where no further calibration scheme is im¬ 
plemented. We see that the errors are reduced for the angle-only 
method when using the ellipticities from IM3SHAPE to measure 
the position angles as compared with the integrated light and mo¬ 
ments based methods. This reduction is due to the correction for the 
PSF when measuring the ellipticities which increases the unifor¬ 
mity of the weighting function and reduces the errors on the angle 
measurements. Hence, the value of w i is lower, and the mul¬ 
tiplicative bias for the zeroth-order estimate is smaller (or equiva¬ 
lently, the /3 WC term as defined in Appendix B is larger) reducing 
the errors on the shear estimates in accordance with equations E 
[4} and ( |B-6[ i. 

The absence of the PSF correction at the measurement stage 
also leads to a larger bias in the zeroth-order shear estimates when 
using the integrated light and the moments based methods. There¬ 
fore, a larger number of iterations are required to remove the bias in 
these methods as shown in equation {32}. The scatter on the shear 
estimates recovered using full ellipticity information is smaller than 
when using only the position angles. This is partly due to form of 
the Fi (|g|) function, as discussed in W14. However, there is also 
an increase in the errors due to the iterative procedure employed 
to remove the biases in the angle-only method; this is quantified in 
equation ( |B-6[ l. 

For all methods using the angle-only estimator, a residual bias 
is expected when using a finite number of iterations (for these 
tests we use a single iteration). However, the limiting factor in this 
method is expected to be the accuracy with which we can simulate 
the galaxies in the branch. 

The results of our analyses of the ege branch are presented 
in Table [I] The error bars on each Q-value were determined by 
assuming that the estimates of m; and a are the true values. We 
used simulations to confirm that the errors on the m, and a es¬ 
timates are approximately Gaussian. We then simulated a suite of 


1 https://github.com/barnabytprowe/great3-public/ 2 http://great3.jb.man.ac.uk/leaderboard/board/ 

wiki/Metric-evaluation-scripts-and-truth-data post_challenge/24 


© 2015 RAS, MNRAS 000.[T]£T3] 

















Angle-only estimators and the GREAT3 simulations 9 




9i 92 




9i 92 

Figure 5. The top panels show the recovered shear estimates for each of the 200 fields of the cgc branch plotted against the input shear 
values and using the integrated light method to measure the position angles. The bottom panels show the recovered shear estimates using 
the moments based method to measure the position angles. In all cases, the black line is the one-to-one line and the red-dashed line is 
the line of best fit for the recovered shear values. 


Method 

Q 

m_|_ (x 10 3 ) 

m x (xlO 3 ) 

c+ (xlO -4 ) 

c x (xlO’ 4 ) 

Int. light 

324 tlH 

-11.26 ±8.25 

-7.34 ±6.22 

-4.03 ± 2.08 

4.34 ± 1.58 

Moments 

403ti° 0 ? 

2.44 ± 7.96 

9.38 ±4.83 

5.32 ± 2.01 

3.32 ± 1.23 

IM3SHAPE: Angle-only 

071 +96 
-*--98 

6.80 ± 3.88 

10.35 ±4.05 

1.89 ±0.98 

1.78 ± 1.03 

IM3SHAPE: Full ellip. 

n7ii! 

29.47 ± 5.18 

27.35 ±4.31 

-11.33 ± 1.31 

0.36 ± 1.10 

IM3SHAPE: Highest entry 

416-50 

0.00 ± 5.00 

-2.09 ±4.14 

-11.0 ± 1.26 

0.44 ± 1.05 

KSB: Highest entry 

12211® 

22.7 ± 7.3 

32.5 ±5.9 

6.19 ± 1.85 

— 1.07 ± 1.51 


Table 1. The results of the analyses we have performed on the GREAT3 cgc branch. The first column shows the Q-values achieved for 
each method. The next two columns are the multiplicative biases estimated for each method. The final two columns show the estimated 
additive biases. 
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Figure 6. The top panels show the shear estimates recovered using the angle-only method for each of the 200 fields of the cgc branch 
plotted against the input shear values and using IM3 SHAPE to measure the position angles of the galaxies. The bottom panels show the 
recovered shear estimates using the full ellipticity information from IM3 SHAPE but with no additional calibration scheme. In all cases, 
the black line is the one-to-one line and the red-dashed line is the line of best fit for the recovered shear values. 


5 x 10 5 uncorrelated Gaussian random variables distributed about 
each of the estimated bias parameters with the dispersions provided 
by the metric evaluation script. These simulated bias parameters 
were used to calculate the 5 x 10° corresponding Q-values using 
the least squares method, and then a histogram of these values was 
constructed. The error bars were calculated as the 68.3% confi¬ 
dence regions of the simulated Q-values distributed about the mean 
recovered value. 

From Table [T] we see that the Q-values for the three angle- 
only analyses are competitive with the highest entry submitted us¬ 
ing IM3SHAPE to within one standard deviation. The Q-values 
achieved in all of the analyses performed with the angle-only ap¬ 
proach are greater than the highest entry submitted using the KSB 
method. Figure [7] displays the values of m; and a calculated us¬ 
ing the GREAT3 metric evaluation script. The plus subscript cor¬ 


responds to the direction parallel with the PSF and the cross sub¬ 
script to the direction perpendicular to the PSF. We see that the 
multiplicative biases of the integrated light and moments based 
angle-only methods lie within two standard deviations of the target 
values identified for future “Stage IV” weak lensing experiments 
aimed at precision dark energy constraints. There appears to be a 
greater residual additive bias when using these methods to mea¬ 
sure the position angles than when using the angles recovered by 
I M3 SHAPE. This is likely to be due to there being no correction 
for the PSF during the angle measurement stage, and hence there 
is a larger additive bias for the zeroth-order shear estimate. Assum¬ 
ing that this is the case, the bias should decrease as the number 
of iterations is increased, as indicated by equation ( [32) , until the 
threshold enforced by the accuracy of the simulations is reached. 
The angle-only method using IM3SHAPE exhibits a reduction in 
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Figure 7. The multiplicative and additive biases for each test on the cgc branch. The grey region indicates the target bias values for 
future precision cosmic shear experiments, which are based on the target values evaluated for the Euclid space mission JMassey et al.| 
|20T3) . 


the multiplicative bias as compared with that of the naive full el- 
lipticity approach, and we also see a reduction in the component of 
the additive bias which aligns with the PSF. 


6 DISCUSSION 

We have demonstrated an algorithm for applying the angle-only 
weak lensing estimator of W14 to realistic galaxy images. We find 
that the performance of this technique is competitive with state- 
of-the-art shape measurement techniques. To measure the position 
angles of the galaxies, we considered three separate approaches: 
the integrated light method, a moments based method, and using 
ellipticities measured by IM3SHAPE. All three methods yield Q- 
values consistent with those of the highest IM3SHAPE submission 
and greater than the highest KSB submission. 

The angle-only estimator requires a measurement of the 
Fi (|g|) function. We have shown that this function can be suc¬ 
cessfully recovered from deep field images using IM3 SHAPE. Con¬ 
straints on the accuracy of the recovered F\ (|g|) function required 
to provide reliable shear estimates are discussed in W14. When 
measuring the ellipticities of the galaxies in the deep field images 
using IM3SHAPE, we necessarily fit for the half-light radii and 
the bulge to total flux ratios of the galaxies. We have introduced 
a method of using this information to construct accurate calibration 
simulations. To do this, we modify the distributions of the fluxes 
and half-light radii measured from the deep fields to match the dis¬ 
tributions in the field being analysed by simply summing over the 
pixel values in that field. This approach could be useful for any 
analysis which requires calibration simulations. It is, as yet, unclear 
as to how general the calibration methods developed in this paper 
are with regards to how the galaxies of a particular survey are se¬ 
lected. A detailed investigation of this issue is left for future work. 
However, provided that one can construct accurate zero shear cal¬ 
ibration simulations, the multiplicative and additive biases present 


in the zeroth-order shear estimates can be corrected for using the 
iterative approach. 

The presence of an anisotropic PSF and pixelization bias po¬ 
sition angle estimates. We have employed a weighting scheme to 
reduce this effect. The weighting function is estimated using cal¬ 
ibration simulations which assume a zero input shear signal. Us¬ 
ing the iterative method employed for this paper, we have argued 
that residual biases in the shear estimates can be reduced below a 
threshold which is determined by the accuracy of the calibration 
simulations. For a perfect suite of calibration simulations, the mul¬ 
tiplicative bias in the zeroth-order shear estimate is essentially due 
to noise in the image, pixelization and the PSF. The magnitude of 
the multiplicative bias in the zeroth-order shear estimate determines 
the rate of convergence for the estimator. Therefore, correcting for 
the PSF prior to the zeroth-order estimate should increase the rate 
of convergence. However, the effect of the PSF on a galaxy’s orien¬ 
tation is model dependent and is, therefore, difficult to correct for 
when using an angle measurement method which is independent of 
the ellipticity. 

The purpose of this paper is to demonstrate the feasibility of 
performing an angle-only shear analysis on realistic galaxy images. 
To demonstrate this, we have focused on the simulated images of 
the cgc branch of the GREAT3 challenge. In future work, we will 
look at how the angle-only method can be applied to fields with 
a variable shear and a variable PSF. We also aim to reduce the de¬ 
pendence of the method on simulations and develop a deeper under¬ 
standing of the level at which the angle-only method can comple¬ 
ment ellipticity based methods and ultimately reduce systematics 
in weak lensing surveys. 
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APPENDIX A: RESIDUAL BIAS IN THE ITERATIVE 
METHOD 


We derive the residual bias in the shear estimates as a function of 
the number of iterations used. 

If we assume the zeroth-order shear estimate can be written as 
equation (|28|>, the first-order iteration gives the shear estimate as 


9i 1] =5 I (0) - (slsim - 9i 0) ) > 
=g ( ° ] - m! gf' ) - c' + Sg^ . 


(A-l) 


We can substitute equation ( |28| ) into equation ( |A-l| l, giving 

9 i 1} -gi-m'i (m'ig + c') + (l - m) Sg^ + SgW. (A-2) 


The second iteration yields 

~(2) -(0) /-(l) / . r (2) 

91 =91 - m i9i - Ci + og\ 

=gt + m? (m-fli + c') + [l - m' (l - to')] 5gf ] 

~ m' i 5gf ) + 5g < j 1) . (A-3) 


From here we see that the error terms propagate linearly through 
the iterations. Hence, as the mean of these terms is zero, there will 
be no residual bias contribution from these terms. If we ignore the 
noise terms, we see that the shear estimate from the n th iteration 
can be written as 

<?i n) = 9i+ (-m'i)" (m'fli + c-) , (A-4) 

and therefore the bias on the n th iteration is 

" .</;) = (-«')" ( m'igi + c') . (A-5) 


APPENDIX B: ERROR ON THE FIRST ITERATION 

We discuss the first-order approximation of the error on the angle- 
only shear estimates with a first-order iteration. For each galaxy 
in the field, we assume that there is an identical galaxy with the 
intrinsic ellipticity rotated by 90°. 
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Let us begin by assuming a first-order approximation of the 
F\ (|g|) function (discussed in W14), such that 

(Isl) « u\g\ , (B-l) 

where u is the first-order coefficient. We write the measured po¬ 
sition angle of a galaxy as equation and assume that errors on 
the position angles are distributed symmetrically about zero. The 
zeroth-order estimate, g^ (with a similar analysis also holding for 
g 2 0> ), can then be written as an average over the galaxy pairs 

JV 

i 2 r 

gi 0) « - ^ Wi(ai + 2^ cos (2a t + 2 
i =1 L 

+ UH ( 2 cti + 7T + 2(5af > ^ cos ^2 on + 7r + 2 

(B-2) 

where there are N/2 galaxy pairs, and where 5ot~A is the error on 
the position angle of the first galaxy in the pair and 5a\ ' is the 
error on the position angle of the corresponding 90° rotated galaxy. 
The weighting function w (at) is normalized and given as 

W (tti) 


W (tti) = 


e r=i «'(«o 


(B-3) 


where the summation is over all galaxies in the field. 

From this form of the estimator, we can write an approximate 
form for the error on the zeroth-order estimate in the limit g —> 0 




*2L( 

2 u 2 \ 


W E cos I 


(2a + 2<k* (1) ) 

2 <Sa (2) ) 


+ w ^ 2) cos (^2a + 7r + 


(B-4) 


where /3 WC is the weighted mean cosine of the error distribution, 


wc — E w ( &i ) cos (2 Sat). 


(B-5) 


We have assumed that the weighting scheme has successfully re¬ 
moved all contributions from an additive bias. From this approach 
we find that cr.( 0 ) ~ <r„(o) = cr a (o). A similar approach can be 
applied to the errors on the shear estimates recovered from the first 
iteration simulations but with the number of galaxies being depen¬ 
dent on the number of galaxies in the simulations; for our analyses 
we simulated ,/V s j m = 10 5 galaxies as discussed in Sectionjd] The 
error on the first iteration of the shear estimator can then be found 
by looking at equation ( | A- 2 | > where the error on the zeroth-order 
shear estimate, g^°\ is Sg[ (as defined in equation (28 ) and the 
error on the estimate from the first-order simulations, < 7 ^, is 
(as defined in equation (|30)). Assuming that the errors on the shear 
estimates are Gaussian (which is expected to be true to first-order 
in the shear due to the central limit theorem), the error on the first- 
order shear estimate is 

<7g fa (i -m') 2 o3 (0) +ct? ( i,, (B- 6 ) 

where a? (0 ) = ^ and = ^dgr^ In the absence of 

an additive bias, the multiplicative bias is effectively due to the /3 C 
term in equation However, when using weighted trigonomet¬ 
ric functions /3c = Pwc, such that m' = /3 WC — 1 . 
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